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SUMMARY 


The  shock  Hugoniot  and  adiabatic  release  curves  of  sedimentary  rocks 
display  a  range  of  different  behaviors  that  affect  the  decay  of  shock  waves 
propagating  away  form  a  confined  source.  Hysteresis  in  the  shock-release  paths  of 
geologic  materials  results  in  irreversible  energy  deposition,  thus  depriving  the 
shock  wave  of  energy  required  to  drive  its  propagation  and  resulting  in  a  more  rapid 
wave  amplitude  decay  and  less  efficient  coupling  to  far  field  seismic  waves  than 
would  be  expected  from  geometrical  effects  alone.  Thus,  data  constraining  this 
behavior  in  rocks  are  essential  for  determining  the  cavity  volumes  for  various 
degrees  of  decoupling  for  different  lithologies. 

We  have  studied  the  shock  compression  and  adiabatic  release  behavior  of 
Coconino  Sandstone  with  an  initial  density  =  2.313  Mg/m3  and  Solenhofen 
Limestone  with  initial  density  =  2.613  Mg/m3,  and  a  limestone  from  Bedford, 
Indiana,  with  initial  density  p^,  =  2.418  Mg/m3.  We  used  a  series  of  shock 
Hugoniot  and  release  adiabat  measurements  on  sandstone  and  limestone  to  develop 
models  for  the  shock  and  release  behavior  of  these  rocks.  The  sandstone  shock 
Hugoniot  curve  can  be  calculated  from  the  equations  of  state  of  quartz,  stishovite, 
and  a  high  coordination  number  liquid.  A  new  fit  of  available  data  for  the  dense 
liquid  phase  gives  p00  =  4.150  Mg/m3,  KSo  =  313.5  GPa,  K'  =  1.8484,  y  = 
1.4(V/V0)0  S,  and  =  2.4  MJ/kg.  In  the  case  of  Coconino  Sandstone,  the 
presence  of  a  few  percent  of  non-quartz  constituents  requires  that  the  STP  densities 
of  the  relevant  phases  be  reduced  by  2.4%  to  fit  the  data.  The  mixed-phase  regions 
of  the  Hugoniot  can  be  modeled  as  a  simple  mixture  of  the  two  relevant  end 
member  phases.  Formation  of  the  higher  pressure  phase  appears  to  be  kinetically 
inhibited,  with  the  fraction  f  of  the  higher  pressure  phase  being  described  by 
f  =  Ae"E,/RT.  For  formation  of  stishovite  from  quartz,  A  =  7.286  and 
Ea(kJ  /  mol)  =  89.36  -  71.97(p0/Poo)  where  poo  and  p0  are  the  actual  and 

theoretical  nonporous  densities  of  the  initial  rock  sample.  For  formation  of  the 
dense  liquid  from  stishovite,  A  =  6.245  x  10^  and  Ea  (kJ/mol)  =  404  +  0.312p00. 

No  formation  of  the  higher  pressure  phase  occurs  until  the  phase  boundary  is 
reached.  As  crude  approximations  we  describe  the  quartz-stishovite  phase 
boundary  by  P  =  12  GPa  and  the  stishovite-dense  liquid  phase  boundary  by 
T  =  4348  -  12.121P  for  T  in  K  and  P  in  GPa. 

Release  adiabats  from  the  single  phase  regions  of  the  sandstone  Hugoniot 
are  modeled  by  the  isentropes  of  the  relevant  phases.  At  ~5  GPa,  stishovite  begins 


iv 


to  form  a  diaplectic  glass.  We  model  the  release  adiabat  below  this  pressure  as  a 
straight  line  in  P-V  space,  from  the  stishovite  release  isentrope  at  5  GPa  to  the  STP 
diaplectic  glass  density  of  2.27  Mg/m3.  The  dense  liquid  undergoes  a  similar 
change  beginning  at  ~5  GPa  with  formation  of  normal  silica  glass  with  an  STP 
density  of  2.204  Mg/m3.  The  quartz  release  adiabat  follows  the  quartz  isentrope  to 
P=0.  As  with  the  Hugoniot,  we  model  the  release  from  the  mixed-phase  regions  of 
the  Hugoniot  as  a  mixture  of  the  relevant  phases,  with  the  proportions  frozen  at  the 
Hugoniot  values. 

The  limestone  Hugoniot  curve  is  complicated  by  several  metastable  phases 
and  multiple-wave  shock  structure  at  low  pressures.  Single  phase  portions  of  the 
Hugoniots  of  Solenhofen  and  Bedford  limestones  are  modeled  as  calcite  I  at  low 
pressures  and  a  dense  liquid  at  high  pressures.  We  have  fit  an  equation  of  state  for 
the  dense  liquid  phase,  with  p0  =  3.377  Mg/m3,  Kso  =  139  GPa,  K  =  3.4,  y  = 

2.5 (V/Vo),  and  E*  =  2.142  MJ/kg.  At  intermediate  pressures,  dense  limestones 
such  as  Solenhofen  seem  to  form  calcite  III,  while  the  more  porous  Bedford 
limestone  apparently  forms  aragonite.  The  mixed  phase  regions  of  the  Hugoniot 
are  modeled  in  a  manner  analogous  to  that  used  for  the  sandstone  Hugoniot.  The 
parameters  describing  the  fraction  f  of  the  higher  pressure  phase  in  the  mixed  phase 
regions  are  A  =  5.848  and  Ea  =  62.7 13-20.651  poo  for  formation  of  calcite  III  or 
aragonite  from  calcite  I  and  A  =  2.919  and  Ea  =  382.453-  130.448Poo  for  formation 
of  the  dense  liquid  from  either  aragonite  or  calcite  III,  for  Ea  in  kJ/mol  and  poo  in 
Mg/m3.  The  phase  boundaries  in  our  model  are  approximated  by  T^  =  164.7P- 
3.6,  Tmh  =  78.00-1.3,  TM.Uq  =  745+48P,  and  Tin.Uq  =  401+25.8P. 

As  with  sandstone,  we  model  the  release  adiabats  with  the  isentropes  of  the 
relevant  phases,  mixing  them  at  the  appropriate  ratios  for  release  from  the  mixed- 
phase  portions  of  the  Hugoniot.  There  is  no  evidence  of  formation  of  a  glass  at 
very  low  pressures,  but  release  of  phases  other  than  calcite  I  is  apparently 
accompanied  by  decomposition  with  evolution  of  C02  gas  at  pressures  below  1 
GPa.  We  have  ignored  this  very  low  pressure  effect  in  our  calculations. 

These  results  allow  us  to  estimate  the  fraction  f  of  shock  internal  energy 
irreversibly  deposited  in  these  rocks  and  thus  no  longer  available  to  drive  the  shock 
wave.  We  find  that  f  is  a  factor  of  ~2  smaller  in  limestones  than  in  sandstone,  so 
that  decoupling  of  the  source  from  far-field  seismic  waves  is  more  efficient  in 
sandstones.  The  value  of  f  may  be  reduced  considerably  in  the  presence  of  an 
elastic  precursor  wave  and  of  intermediate  waves.  Experimental  evidence  indicates 
that  shock  wave  decay  still  occurs  even  in  the  supposedly  purely  elastic  shock 


regime,  suggesting  that  effects  such  as  fracture  formation  are  very  important  to  the 
overall  picture  of  explosive  energy  coupling. 
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I.  INTRODUCTION 


The  shock  Hugoniot  and  adiabatic  release  curves  of  sedimentary  rocks 
display  a  range  of  different  behaviors  that  affect  the  decay  of  shock  waves 
propagating  away  form  a  confined  source.  Many  of  the  minerals  making  up 
sandstones,  shales,  and  limestones  undergo  phase  transitions  when  shocked  to 
pressures  of  interest  to  studies  of  coupling  of  energy  from  explosive  sources  into 
far  field  seismic  waves.  Both  sandstones  and  limestones  have  been  observed  to 
exhibit  elastic  precursors  and  multiple  wave  behavior  at  shock  velocities  up  to  3.7 
km/s  and  5.7  km/s,  respectively  [Ahrens  and  Gregson,  1964].  Hysteresis  in  the 
shock-release  paths  of  materials  results  in  irreversible  energy  deposition,  thus 
depriving  the  shock  wave  of  energy  required  to  drive  its  propagation  and  resulting 
in  a  more  rapid  wave  amplitude  decay  and  less  efficient  coupling  to  far  field  seismic 
waves  than  would  be  expected  from  geometrical  effects  alone.  Thus,  data 
constraining  this  behavior  in  rocks  are  essential  for  determining  the  cavity  volumes 
for  various  degrees  of  decoupling  for  different  lithologies. 

We  have  studied  the  shock  compression  and  adiabatic  release  behavior  of 
Coconino  Sandstone  with  a  mean  initial  density  of  p^,  =  2.313  Mg/m3  and 
Solenhofen  Limestone  with  a  mean  initial  density  of  Pqq  =  2.613  Mg/m3,  and  a 
limestone  from  Bedford,  Indiana,  with  a  mean  initial  density  of  p00  =  2.418  Mg/m3. 
The  sandstone  is  primarily  quartz,  with  ~3%  feldspar.  The  limestones  are 
essentially  pure  calcite  with  very  minor  amounts  of  other  constituents. 
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n.  SHOCK  WAVE  EQUATION  OF  STATE  OF  SANDSTONE  AND 

LIMESTONE 

A.  Experimental  Techniques. 

The  experimental  techniques  used  for  this  study  are  discussed  in  detail 
elsewhere  [Ahrens,  1987].  A  sample  of  the  rock  being  studied  is  placed  on  a  metal 
driver  plate  (Figure  1),  with  standard  buffers  and/or  inclined  mirrors  placed  on  the 
free  (rear)  surface  of  the  sample.  Shock  arrival  mirrors  are  placed  against  the  rear 
surfaces  of  the  driver  plate,  sample,  and  buffers.  During  the  experiment,  the  driver 
plate  is  impacted  by  a  gun-launched  projectile  containing  a  metal  impactor  (flyer 
plate).  The  impact  generates  a  shock  wave  which  propagates  through  the  driver 
plate,  sample,  and  buffer  materials.  Arrival  of  the  shock  wave  at  the  free  surface  of 
each  component  of  the  target  causes  disruption  of  the  specular  reflection  from  the 
shock  arrival  mirrors  as  shown  in  Figure  lb.  One  spatial  dimension  of  the  target, 
including  all  shock  arrival  mirrors,  is  imaged  as  a  function  of  time  by  a  streak 
camera,  allowing  the  shock  wave  transit  times  and  thus  the  shock  wave  velocities 
through  target  components  of  known  thickness  to  be  determined.  The  internal 
energy  E,  pressure  P,  density  p,  and  particle  velocity  up  of  the  shocked  sample  and 
buffer  are  calculated  from  the  measured  projectile  velocity,  sample  initial  density, 
and  shock  wave  velocity  using  the  Rankine-Hugoniot  shock  equations  for 
conservation  of  mass,  momentum,  and  energy  and  by  the  requirement  that  pressure 
and  particle  velocity  must  be  continuous  across  material  interfaces. 

The  buffers  and  inclined  mirrors  are  used  to  measure  points  on  the  release 
adiabats  of  the  samples  by  the  requirement  that  the  P-Up  state  in  the  released  state  of 
the  sample  be  identical  with  that  in  the  shocked  buffer,  or  the  free  surface  in  the 
case  of  the  inclined  mirror.  Table  1  give  the  properties  of  the  materials  used  for 
flyer  and  driver  plates  and  buffers  in  this  study. 

The  shock  Hugoniot  of  a  material,  which  is  the  locus  of  final  states  achieved 
by  shock  compression,  is  conveniently  described  as  a  polynomial  representation  of 
the  shock  wave  velocity,  Us  in  terms  of  the  particle  velocity,  ty 

Us=C0+sup+s'up2  (1) 
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Shock  Arrival 
at  Driver-sample 
Interface 


Shock  Arrival  at 
Sample  Surface 


Time 


Shock  Arrival 
at  Buffer  Surface 


Steak  Camera  Slit  Image 


Figure  1.  Diagram  of  a  typical  shock  Hugoniot  experiment  with  buffers  for  release 
adiabat  measurements,  (a)  Projectile  in  flight  prior  to  impact,  (b)  Static  (left)  and 
dynamic  image  recorded  by  streak  camera. 
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Table  1.  Shock  Hugoniot  parameters  for  standard  materials  used  in  this  study. 

Material  rO  (Mg/m3)  CO  (km/s)  s  s’ (s/km)  Source 


All  100 

2.707 

5.386 

1.339 

A12024 

2.784 

5.330 

1.34 

W 

19.224 

4.029 

1.237 

Ta 

16.650 

3.293 

1.308 

Lexan 

1.193 

2.449 

1.498 

4.817 

0.589 

1.907 

1.431 

Polystyrene 

0.0474 

0.362 

0.9165 

0.04982 

(1)  Marsh  [1980] 

(2)  Mitchell  and  Nellis  [1981] 

(3)  Fit  to  data  from  Marsh  [1980]  for  up  <  2.605km  /  S. 

(4)  Fit  to  data  from  Marsh  [1980]  for  2.605km  /  S  <  up  <  3.456km  /  S. 

(5)  Fit  to  data  from  Marsh  [1980]  for  Up  >  3.456km  /  S. 

(6)  Calculated  from  p00  using  expression  from  Anderson  and  Ahrens  [1993] 
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Usually,  the  polynomial  is  truncated  at  first  or  second  order.  The  pressure  P, 
density  p,  and  internal  energy  E  in  the  shocked  state  are  given  by 

Ph  =  PooUsUp 


(2) 


Ph  —  Poo 


Us 

Ut-up 


(3) 


Here,  the  subscripts  0  and  H  refer  to  the  unshocked  and  shocked  (Hugoniot)  states, 
respectively,  and  the  subscript  00  on  p  refers  to  the  initial  density ,  including  any 
porosity  in  the  sample.  We  arbitrarily  set  E  =  0  as  the  initial  condition. 


B.  Sandstone 
1.  Experimental  Data 

Table  2  presents  the  experimental  data  for  the  Coconino  Sandstone  used  in 
this  study,  which  had  a  mean  initial  bulk  density  of  2.313  Mg/m3.  The  densities  of 
the  release  states  are  upper  limits,  obtained  by  approximating  the  release  paths  as 
straight  line  segments  in  the  P-V  plane  [Lyzenga  and  Ahrens,  1978]. 

Figure  2  shows  the  present  data,  along  with  an  empirical  fit  to  the  Us-up 
projection  of  the  Hugoniot  The  Us-up  projection  of  the  Hugoniot  curve  can  be 
divided  into  two  straight-line  segments  with  weighted  fits  of  C0  =  3.210+.026  km/s 
and  s  =  0.978±.01 1  for  up  <  2.932  km/s  and  C0  =  0.973+.108  km/s  and  s  = 
1.741+036  for  up  >  2.939  km/s.  Also  shown  is  the  P-V  projection  of  the  curve 
described  by  the  empirical  Us-Up  fit,  along  with  the  present  data  and  data  from  Van 

Thiel  [1977]  for  comparison. 

Sandstones  generally  show  an  elastic  precursor  for  very  low  shock 
pressures.  Based  on  the  data  from  Ahrens  and  Gregson  [1964],  the  velocity  Ud  of 
the  elastic  precursor  is  only  very  weakly  dependent  on  the  particle  velocity,  with 
even  the  sign  of  the  variation  being  uncertain.  As  a  first  order  approximation,  the 
elastic  precursor  velocity  depends  only  on  the  initial  density  of  the  rock .  Using  the 
data  from  Ahrens  and  Gregson  [1964],  which  cover  polycrystalline  quartz  rocks 
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Table  2.  Shock  wave  and  release  adiabat  data  for  Coconino  Sandstone. 
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Figure  2.  Shock  Hugoniot  and  release  adiabat  data  and  Hugomot  based  on 
empirical  Us-up  fits  for  Coconino  Sandstone,  (a)  Us-up  plane,  (b)  P-V  plane. 
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with  poo  in  the  range  from  1.961  Mg/m3  to  2.650  Mg/m3,  we  describe  the  elastic 
precursor  velocity  Ue)  by 

Uel=4.53p00  -6.122  (5) 

where  Uel  is  given  in  km/s  and  Pqq  is  in  Mg/m3. 

2.  Equation  of  State. 

Although  the  experimental  Us-Up  relation  gives  us  a  good  empirical 
description  of  the  Hugoniot  of  the  particular  rock  we  studied,  we  wish  to  develop  a 
model  for  shock  and  release  behavior  that  is  generally  applicable  to  quartz  rocks  in 
general.  Such  a  model  could  be  easily  modified  to  account  for  minor  non-quartz 
constituents. 

Since  the  predominant  constituent  is  quartz,  the  phases  appearing  on  the 
quartz  shock  Hugoniot  curve  are  expected  to  also  appear  on  the  sandstone  Hugomot 
curve.  In  particular,  these  are  quartz,  stishovite,  and  dense  Si(>2  liquid.  Above  the 
Hugoniot  elastic  limit  (HEL),  pressure  at  a  given  specific  volume  VH  =  l/pH  on  the 
Hugoniot  can  be  calculated  by 

p  _  Es  +  Ett  —  PsVh/Yh  (6) 

H  (Voo-Vh)/2-Vh/yh 

where  PH  is  the  shock  state  (Hugoniot)  pressure,  Ps  is  the  pressure  of  the  material 
on  the  principal  isentrope  at  V  =  VH,  VH  is  the  shock  state  specific  volume,  is 
the  initial  specific  volume  of  the  material,  Ett  is  the  phase  transition  energy  at  STP, 
Es  is  the  internal  energy  gained  in  compression  of  the  material  along  the  principal 
isentrope  to  specific  volume  VH,  and  Yh  =  Griineisen  parameter  of  the  material  at 
volume  VH.  This  expression  equates  the  energy  at  pressure  PH  and  specific  volume 
VH  obtained  by  following  the  thermodynamic  path  of  the  shock  compression  (a 
straight  line  in  P-V  space)  with  that  obtained  via  phase  change  (if  needed)  at  room 
temperature  and  1  bar  pressure,  isentropic  compression  to  volume  VH,  and 
isochoric  heating  to  pressure  PH.  We  assume  that  the  principal  isentrope  is 
described  by  the  third  order  Eulerian  finite  strain  Birch-Mumaghan  equation: 

ps  =  |kso(x,-xsXi+5-^!)  (7) 
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5  =  |(4-K') 


(10) 


where  is  the  value  of  the  isentropic  bulk  modulus  Ks  at  STP  and  K’  in  the 
isentropic  pressure  derivative  of  Ks  at  STP.  We  also  assume  that  the  Griineisen 

parameter  y  obeys 

Y  =  T,(V/V0)“  * 


For  the  properties  of  quartz,  we  use  the  equation  of  state  (EOS)  parameters 
obtained  by  averaging  the  values  given  by  Sumino  and  Anderson  [1984].  For 
simplicity,  we  assume  that  Cv=3k/atom.  For  stishovite,  we  use  the  EOS 
parameters  tabulated  by  Lyzenga  et  al.  [1983],  but  with  a  slightly  different  value  of 
Ks0,  which  seems  to  better  describe  the  shock  wave  compression  data  for  quartz 
given  by  Wackerle  [1962]  and  presented  graphically  by  Chhabildas  [1986]  and 
Lyzenga  et  al.  [1983].  The  shock  temperature  data  of  Lyzenga  et  al.  [1983]  suggest 
that  Si02  is  liquid  above  117  GPa  and  70  GPa  on  the  quartz  and  fused  quartz  shock 

Hugoniots,  respectively.  Since  our  data  are  for  samples  which  are  only  slightly 
more  dense  than  fused  quartz  and  the  Coconino  Sandstone  data  tabulated  by  Van 
Thiel  [1977]  are  for  densities  lower  than  that  of  fused  quartz,  we  can  expect  our 
highest  pressure  data  to  be  for  shock  states  in  the  liquid  phase.  The  densities  of  the 
shock  states  for  quartz  and  fused  quartz  in  the  relevant  pressure  ranges  are 
incompatible  with  the  known  properties  of  liquid  Si02  at  low  pressure  and  strongly 

suggest  that  the  liquid  phase  on  the  Hugoniot  is  characterized  by  Si4+  in  6-fold 
coordination  with  O^-,  as  it  is  in  stishovite.  This  is  a  reasonable  conclusion,  since 
stishovite  is  the  solid  which  would  undergo  equilibrium  melting  to  give  this  liquid. 
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Lyzenga  et  al.  [1983]  give  a  set  of  EOS  parameters  for  the  dense  liquid  phase,  but 
we  cannot  reconcile  the  resulting  shock  compression  curves  with  all  the  existing 
Hugoniot  data  for  quartz  and  fused  quartz.  We  have  refit  the  values  of  Ks0  and  K' 
to  data  from  Lyzenga  et  al.  [1983]  for  quartz,  Trunin  et  al.  [1971a]  for  porous 
quartz,  Trunin  et  al.  [1971b]  for  quartz  and  quartzite  and,  Lyzenga  et  al.  [1983]  and 
Marsh  [1980]  for  fused  quartz.  We  retain  the  empirical  expression  for  Cv  given  by 
Lyzenga  et  al.  [1983].  The  resulting  fit  is  shown  in  Figure  3.  The  EOS  parameters 
used  for  the  three  phases  of  interest  to  this  study  are  listed  in  Table  3,  along  with 
slight  variations  of  these  values  which  give  a  better  fit  to  the  data  for  Coconino 
Sandstone.  The  primary  change  made  to  bring  the  model  in  line  with  the  data  was 
to  vary  the  values  of  po  for  the  different  phases.  This  is  presumed  to  reflect  the 
effects  of  non-quartz  constituents  of  the  sandstone. 

Above  the  HEL,  the  phase  on  the  sandstone  and  quartz  Hugoniots 
apparently  remains  quartz  until  sufficient  energy  is  available  to  initiate  transition  to 
stishovite.  Although  part  of  the  Hugoniot  is  in  the  coesite  stability  field,  the 
available  data  show  no  evidence  of  coesite  formation  in  the  shock  state.  A  wide 
mixed-phase  (MP)  region  exists  on  the  Hugoniot  between  the  onset  and  completion 
of  stishovite  formation.  The  lower  and  upper  pressure  boundaries  to  this  region 
depend  on  the  value  of  p00.  The  lower  boundary  on  the  MP  regime  occurs 
considerably  above  the  pressure  at  which  the  Hugoniot  enters  the  stishovite  stability 
field,  but  decreases  as  p^  decreases,  suggesting  that  the  transition  is  kinetically 
inhibited  (see,  e.g.,  Podurets  and  Trunin  [1987]).  We  can  model  the  fraction  f  of 
the  higher-pressure  phase  by  assuming  an  Arrhenius  rate  constant  for  the  transition 
and  that  the  timescale  for  the  transition  is  relatively  constant  (i.e.,  that  the  thickness 
of  the  shock  front  is  proportional  to  Ut): 


f  =  Ae-E*/RT 


(12) 


where  the  preexponential  factor  A  and  the  activation  energy  Ea  are  determined 
empirically  from  the  experimental  data.  The  assumption  here  is  that  no  transition 
occurs  until  the  phase  boundary  is  crossed,  and  then  is  controlled  by  (12)  after  the 
phase  boundary  is  crossed.  The  issue  of  kinetics  raises  the  possibility  that 
specimen  size  could  affect  the  behavior  of  the  Hugoniot.  However,  this  was 
experimentally  tested  by  Podurets  et  al.  [1976],  who  concluded  that  specimen  size 
does  not  affect  shock  propagation  or  the  position  of  the  P- V  Hugoniot  curve.  For 
the  present  model,  we  assume  that  the  quartz-stishovite  transition  occurs  at  a  set 
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Figure  3.  Equation  of  state  fits  to  shock  Hugoniot  data  for  quartz  and  fused  quartz  in  the  high  coordination  number  liquid 
stability  region. 


Table  3.  Equation  of  state  properties  of  phases  for  pure  SiC>2  and  Coconino  Sandstone. 
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pressure  of  12  GPa.  To  simplify  the  calculation  from  a  practical  standpoint,  we 
also  assume  that  the  relevant  temperature  T  is  that  for  the  lower  pressure  phase  at 
the  correct  shock  pressure.  From  the  present  data  and  data  presented  by  Chhabildas 
[1986]  for  quartz  and  Van  Thiel  [1977]  for  Coconino  Sandstone  with  =  1.98 
Mg/m3,  we  find  that  the  preexponential  factor  has  a  constant  value  of  A  =  7 .286, 
while  the  activation  energy  obeys 

Ea(kJ  /  mol)  =  89.36  -  7L97(p0/p00)  (13) 

where  p0  is  the  value  for  quartz. 

The  MP  region  is  modeled  by  a  volume  weighted  average  of  the  Hugoniot 
curves  of  the  low  and  high  pressure  phases  in  the  P-V  plane: 

VMP=fVHpp  +  (l-f)VLpp  <14> 

where  the  subscripts  LPP,  MP,  and  HPP  refer  to  the  low  pressure  phase,  mixed 
phase,  and  high  pressure  phase,  respectively.  This  is  not  thermodynamically 
rigorous,  but  is  convenient  for  the  calculations  and  is  a  very  good  approximation, 
considering  that  thermal  effects  on  V  are  fairly  small  at  high  pressures. 

The  transition  region  between  stishovite  and  the  high-density  liquid  can  be 
treated  in  the  same  way  as  that  between  quartz  and  stishovite.  The  metastable 
Hugoniot  curves  in  this  region  are  very  similar  for  the  cases  where  data  exist  (i.e., 
on  the  quartz  and  fused  quartz  Hugoniot  curves)  so  that  the  transition  region  is  not 
readily  apparent  in  the  P-V  data.  However,  the  transition  region  is  very  apparent  in 
the  P-T  data  of  Lyzenga  et  al.  [1983].  We  use  a  straight  line  fit  to  our  interpreted 
phase  boundary  positions  from  Lyzenga  et  al.  [1983]  as  the  crudest  representation 
of  the  stishovite-liquid  phase  boundary: 

T  =  4348-12.121P  (15) 

for  T  in  K  and  P  in  GPa.  Their  data  show  that,  similar  to  the  quartz-stishovite 
transition,  the  stishovite-liquid  transition  does  not  begin  until  the  Hugoniot  curve  is 
well  into  the  liquid  stability  region.  The  transition  then  occurs  over  a  fairly  narrow 
pressure  range.  The  details  of  the  transition  region  cannot  be  determined  with 
confidence  from  the  existing  data,  so  we  have  chosen  to  express  the  fraction  of 
liquid  by  (12),  with  A  =  6.245  x  109  and  an  activation  energy  given  by 
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Ea(kJ  /  mol)  =  404  +  0.312p00 


(16) 


for  p00  in  kg/nA 

Figure  4  shows  the  Hugoniot  curves  we  calculate  for  quartz  and  Coconino 
sandstone  using  the  EOS  parameters  in  Table  3  and  the  kinetics  model  discussed 
above. 


3.  Release  Adiabats. 

The  release  adiabats  of  each  of  the  Hugoniot  phases  can  be  calculated  from 
the  EOS  parameters  listed  in  Table  3,  under  the  assumption  that  the  constant  volume 
specific  heat,  Cv,  does  not  vary  significantly,  which  is  a  reasonable  approximation. 
For  a  given  density  pr  on  the  isentrope  anchored  to  a  shock  state  at  density  pH  and 
shock  pressure  PH,  the  pressure  Pr  on  the  release  isentrope  is 

P,  =  Ps,  +  (P„  -  PsXvn/V,r'  expjjj"  (7/V)dv]  (17) 


where  PSr  is  the  pressure  in  the  principal  isentrope  corresponding  to  V  =  Vr  in  (7)  - 
(9).  The  internal  energy  difference  between  this  and  the  Hugoniot  state  (i.e.,  the 
energy  converted  from  internal  to  kinetic  energy  by  doing  work  to  expand  the 
material  during  the  release  process)  is 


EH"Er 


(18) 


where  Eh  is  given  by  (4)  The  internal  energy  not  converted  into  kinetic  energy  upon 
release  is  deposited  as  heat.  This  irreversible  energy  deposition  causes  the  shock 
wave-release  wave  couple  to  decay  as  energy  required  to  drive  the  wave 
propagation  is  sapped  away. 

The  particle  velocity  of  material  released  from  the  shocked  state  is  related  to 
the  P-V  release  path  via  the  Riemann  integral  formula  [Rice  et  al.,  1958]: 


(19) 
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Figure  4.  Model  Hugoniots  and  release  adiabats  for  (a)  quartz  and  (b)  Coconino 
Sandstone. 
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where  the  subscript  rel  refers  to  the  isentropic  release  path.  An  additional 
contribution  to  Up,-  can  come  from  the  energies  of  transition  associated  with  phase 
changes  during  release.  We  have  assumed  that  the  total  internal  energy  is  constant 
with  latent  heat  affecting  sensible  heat  rather  than  macroscopic  kinetic  energy.  The 
simplest  model  assumes  that  the  relative  fractions  of  the  low  and  high  pressure 
phases  are  kinetically  frozen  at  their  Hugoniot  state  values.  The  available  evidence 
from  recovered  shocked  quartz  and  from  the  present  data  is  that  the  stishovite 
transforms  to  a  diaplectic  glass  with  an  STP  density  of  2.27  Mg/m3  and  that  dense 
amorphous  SiC>2  transforms  to  the  more  familiar  low-pressure  glass  (p  =  2.204 

Mg/m3)  or  liquid  at  low  pressures.  Based  on  the  available  data,  we  assume  this 
happens  at  a  critical  pressure  Pc  of  ~5.1  GPa  for  both  phases.  In  both  cases,  we 
have  assumed  that  the  path  below  Pc  is  a  straight  line  in  P-V  space  from  the  volume 
at  Pc  to  the  STP  volume  of  fused  quartz.  This  ignores  the  effects  of  thermal 
expansion,  but  the  error  introduced  is  insignificant.  As  with  the  Hugoniot,  we 
model  the  release  points  on  the  mixed-phase  release  paths  by  assuming  a  mass- 
weighted  average  of  the  volumes  of  the  low  and  high  pressure  phases  via  (14). 
Figure  4  presents  release  adiabats  calculated  with  this  model. 


C.  Limestone. 

1.  Experimental  Data. 

Table  4  presents  the  data  obtained  for  Solenhofen  and  Bedford  limestones  in 
this  study.  Figures  5  and  6  present  these  data  along  with  data  from  Ahrens  and 
Gregson  [1964],  Van  Thiel  [1977],  and  Tyburczy  and  Ahrens  [1986]. 

The  low  pressure  behavior  of  CaC03  is  very  complicated.  At  least  four 

phases  are  known  which  seem  to  be  only  slightly  different  in  internal  energy.  In 
addition  to  calcite  I  and  aragonite,  at  least  two  phases,  calcite  II  and  calcite  HI,  are 
observed  in  static  compression,  with  the  calcite  1-13  and  D-III  transitions  occurring 
at  -1.5  GPa  and  2.2  GPa,  respectively,  at  298  K  [Fiquet  et  al.,  in  press].  Static 
compression  data  of  Adadurov  et  al.  [1960]  also  show  a  very  slight  volume 
discontinuity  at  —1.16  GPa,  although  this  could  be  due  to  some  experimental  effect 
Complicating  this  issue  is  that,  with  the  exception  of  calcite  I,  different  investigators 
disagree  markedly  in  their  determinations  of  the  properties  of  these  phases. 

Low  pressure  shock  wave  data  for  calcite  and  calcite  rocks,  including  the 
data  obtained  in  the  present  study,  show  not  only  an  elastic  precursor  wave  at  low 
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Table  4.  Shock  wave  and  release  adiabat  data  for  Solenhofen  and  Bedford  limestones. 
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Bedford  2.4226  A12024  1.5427  0.7472c  4.7672c  8.63c  2.8729c 

(0.0025)  (0.0072)  (0.0105)  (0.1238)  (0.26)  (0.0160) 

0.9372  3.9975  10.40  3.0512  2.81088  0  1.5033 

(0.0113)  (0.0413)  (0.20)  (0.0025)  (0.0576)  (0.0378) 


Table  4.  (continued) 
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Figure  5.  Shock  Hugoniot  and  release  adiabat  data  and  Hugoniot  based  on 
empirical  Us-up  fits  for  Solenhofen  Limestone,  (a)  Us-up  plane,  (b)  P-V  plane. 
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Particle  Velocity  (km/s) 


Figure  6.  Shock  Hugoniot  and  release  adiabat  data  and  Hugoniot  based  on 
empirical  Us-up  fits  for  limestone  from  Bedford,  Illinois,  (a)  Us-up  plane,  (b)  P-V 
plane. 
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shock  stresses,  but  successive  formation  of  at  least  two  and  perhaps  three  high- 
pressure  phases,  detected  by  the  sequential  arrival  of  multiple  shock  fronts. 

Although  the  data  tend  to  be  somewhat  scattered,  the  stress  levels  of  these  different 
waves  in  the  multiple-wave  experiments  are  similar  to  those  characterizing  calcite  II 
and  III  formation  in  static  compression  experiments. 

As  with  sandstone,  the  velocity  of  the  elastic  precursor  wave  is  not  strongly 
correlated  with  particle  velocity  and  is,  to  first  order,  a  function  only  of  the  initial 
density  of  the  limestone.  We  have  combined  the  data  of  Ahrens  and  Gregson 
[1964]  for  the  first  wave  velocities  in  polycrystalline  calcite  rocks  with  those  from 
our  experiments  to  obtain  a  data  set  for  rocks  with  poo  in  the  range  from  2.327 
Mg/m3  to  2.619  Mg/m3.  A  fit  of  Ud  as  a  function  of  p00  gives 

Uel  =-80.256  + 63. 217p00  -11. 617pM2  (2°) 

where  Uel  is  in  km/s  and  p00  is  in  Mg/m3. 

2.  Equation  of  State. 

Because  of  the  confused  nature  of  the  CaCC>3  system  at  low  pressures,  we 
can  provide  only  a  very  approximate  treatment  in  that  region.  One  very  important 
complication  is  that  the  final  shock  states  of  calcite  rocks  below  ~25  GPa  cannot  be 
matched  by  any  published  equation  of  state  for  any  single  phase  of  CaC03-  We 

have  chosen  to  model  the  lowest  pressure  portions  of  the  Hugoniots  of  low- 
porosity  calcite  rocks  as  a  mixture  of  calcite  I  and  calcite  HI,  where  we  identify 
calcite  m  as  the  phase  observed  by  both  Fiquet  et  al.  [1993]  and  Vizgirda  and 
Ahrens  [1982].  In  the  case  of  Bedford  limestone,  (porosity  ~  12%),  aragonite  is 
used  instead  of  calcite  m,  as  explained  below. 

The  identity  of  the  phase  comprising  the  final  shock  state  just  above  this 
mixed  phase  region  is  uncertain.  The  shock  compression  data  for  single  crystal 
calcite  from  Ahrens  and  Gregson  [1964]  seem  to  indicate  that  the  final  phase  is  not 
calcite  HI,  but  a  higher  pressure  phase  called  calcite  IV  by  Tyburczy  and  Ahrens 
[1986].  However,  Fiquet  et  al.  [1993]  present  new  X-ray  diffraction  data  which 
indicate  that  calcite  m  has  the  same  properties  as  calcite  IV  given  by  Vizgirda  and 
Ahrens  [1982].  The  previous  EOS  parameters  for  calcite  III  are  called  into  question 
because  the  phase  coexistence  regions  noted  by  Fiquet  et  al.  [1993]  can  introduce 
errors  into  non-X-ray  volume  determinations  used  in  the  earlier  studies.  Fiquet  et 
al.  [1993]  do  note  the  appearance  of  new  X-ray  lines  and  a  possible  calcite  IV  phase 
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at  about  the  same  pressure  as  found  by  Ahrens  and  Gregson  [1964],  but  the 
properties  of  this  phase  under  shock  compression  seem  to  be  consistent  with  calcite 
HI,  so  the  situation  is  at  present  confused.  A  further  complication  is  that  the 
Bedford  Limestone  Hugoniot  above  this  point  is  more  consistent  with  the  properties 
of  aragonite  than  calcite  HI.  The  issue  will  likely  require  further  careful  static  X-ray 
diffraction  experiments  at  elevated  temperature  and  pressure  to  map  out  the  phase 
diagram  before  it  can  be  resolved. 

Based  on  the  available  data  and  the  EOS  parameters  presented  by  Fiquet  et 
al.  [1993]  and  Vizgirda  and  Ahrens  [1982],  we  assume  shock  states  in  the  25-40 
GPa  range  for  calcite  to  be  a  single  phase  which  we  identify  as  calcite  HI.  Based  on 
the  published  data  and  previous  determinations  of  EOS  parameters,  we  propose  that 
calcite  HI  is  described  by  p^,  =  3.100  Mg/m3,  Kso  =  96.5  GPa,  K'  =  4.0,  y0  = 
1.183,  n  =  0.6,  and  E*  =  0.02  MJ/kg.  These  parameters  also  describe  the 
equivalent  portions  of  the  aragonite  and  relatively  low  porosity  limestone  Hugoniots 
reasonably  well.  For  limestones  with  porosities  as  much  as  15%,  the 
corresponding  region  is  better  fit  by  aragonite. 

Above  the  calcite  HI  and  aragonite  single  phase  regions,  the  Hugoniot  state 
appears  to  change  to  yet  another  phase.  We  fit  the  shock  compression  data  for  P  > 
100  GPa  for  calcite  and  limestones  to  obtain  p0  =  3.377  Mg/m3,  Kso  =  137  GPa, 

K*  =  3.4,  y0  =  2.5,  n  =  1.0,  and  Etf  =  2.142  MJ/kg.  We  presume  that  this  phase 
probably  is  a  dense  melt,  given  the  large  STP  transition  energy  from  calcite  I. 

Table  5  summarizes  the  EOS  parameters  we  have  adopted  for  each  phase 
considered  in  this  study.  We  assumed  in  all  cases  that  Cv  =  3k/atom.  Calcite  II  is 
not  actually  used  in  our  study,  but  is  presented  for  completeness.  Another 
complication  with  calcite  rocks  is  that  no  combination  of  the  phases  described  in 
Table  5  will  give  the  correct  Hugoniots  for  highly  porous  rocks  such  as  chalk  (see, 
e.g.,  Kalashnikov  et  al.  [1973],  and  Tyburczy  and  Ahrens  [1986]).  Because  of  the 
large  internal  energies  induced  by  shock  waves  in  porous  materials,  it  is  probable 
that  the  phases  encountered  on  the  Hugoniots  of  these  materials  are  not  the  same  as 
those  found  on  the  Hugoniots  of  the  denser  rocks  studied  here.  We  attempted  to  fit 
separate  equations  of  state  for  the  data  for  chalk,  but  still  could  not  get  consistent 
results.  We  suggest  that  the  Hugoniot  states  of  very  porous  calcite  rocks  may  be 
affected  by  decomposition  of  the  CaCC>3.  Such  behavior  will  exacerbate  the 
difficulties  of  developing  accurate  models  of  the  behavior  of  porous  CaCC>3  rocks. 
At  present,  we  advise  that  the  models  developed  in  this  study  be  restricted  to 
limestones  with  <20%  porosity. 
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Table  5.  Equation  of  state  properties  of  CaCC>3  phases. 


Property 

Calcite  I 

Calcite  n 

Aragonite 

Calcite  HI 

Dense  Liquid 

p0  (Mg/m3) 

2712 

2522 

2931 

3100 

3377 

KS0(GPa) 

74.405 

13 

66.73 

96.5 

139 

K' 

5.42 

4 

4 

4 

3.4 

yo 

0.554 

1.5 

1.80 

1.183 

2.5 

n 

1 

1 

1 

0.5 

1 

Etr  (MJ/kg) 

0 

0 

-0.0006 

0.02 

2.142 

Cv  (J/kg/K) 

1246 

1246 

1246 

1246 

1246 

Source 

(1) 

(2) 

(3) 

(4) 

(5) 

(1)  Sumino  and  Anderson  [1984] 

(2)  Fit  to  data  from  Fiquet  et  al.  [1993] 

(3)  Calculated  from  Salje  and  Viswanathan  [1976]  and  Robie  et  al.  [1978] 

(4)  Fiquet  et  al.  [1993]  and  Vizgirda  and  Ahrens  [1982] 

(5)  This  study 
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As  with  the  SiC>2  system,  we  model  the  mixed-phase  regions  of  the 
Hugoniot  by  assuming  kinetically-inhibited  transformations  on  the  equilibrium 
phase  boundary  is  crossed.  The  primary  reason  for  taking  this  assumption  is  the 
fact  that  formation  of  the  high  pressure  phase  identified  as  calcite  in  by  Fiquet  et  al. 
[1993]  above  2.2  GPa  is  not  complete  until  ~25  GPa  for  single-crystal  calcite  under 
shock  compression.  Using  the  EOS  parameters  in  Table  5  to  calculate  the  Hugoniot 
temperatures  corresponding  to  the  P-V  data,  we  fit  Arrhenius  parameters  for  (12)  to 
get  A  =  5.848  and 

E.  =  62.713 -20. 65  lpoo  <21) 

for  pOT  in  Mg/m3  and  Ea  in  kJ/mol  for  the  conversion  of  calcite  I  to  calcite  HI  or 
aragonite.  For  the  calcite  ID  -  liquid  and  aragonite  -  liquid  mixed  phase  regions,  we 
get  A  =  2.919  and 

Ea  =  382.453  -  130.448poo  (22) 


for  poo  in  Mg/m3  and  Ea  in  kJ/mol.  For  the  calcite  I-aragonite  and  calcite  I-calcite  m 
phase  boundaries,  we  take  the  overall  trends  of  the  data  presented  by  Jamieson 


[1957]: 

Ti-Ar=164-7p-3-6 

(23) 

TI_m=78.0P-1.3 

(24) 

for  T  in  K  and  P  in  GPa.  For  the  provisional  calcite  IE-liquid  and  aragonite-liquid 
phase  boundaries,  we  get 

T*_u,  =  745  +  48P 

Tm_uq  =  401  +  25. 8P 

for  T  in  K  and  P  in  GPa.  These  phase  boundaries  are  extremely  uncertain,  being 
based  solely  on  the  behavior  of  the  P-V  Hugoniots  and  the  temperatures  calculated 
from  the  EOS  parameters  given  in  Table  5.  Figure  7  presents  the  Hugoniots  of 
Solenhofen  and  Bedford  limestones  calculated  with  this  model. 


(25) 

(26) 
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Figure  7.  Model  Hugoniots  and  release  adiabats  for  (a)  Solenhofen  and  (b) 
Bedford  limestones. 
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3.  Release  behavior. 

Measured  points  on  release  adiabats  from  the  calcite  in  portion  of  the 
Hugoniot  seem  to  require  reversion  to  calcite  I  with  energy  approximately  equal  to 
the  STP  transition  energy  being  converted  to  kinetic  energy.  One  problem  with  this 
interpretation,  however,  is  that  the  phase  reversion  takes  place  at  a  pressure  higher 
than  would  be  expected  based  solely  on  the  CaC03  phase  transition  pressures.  We 
may  be  detecting  partial  decomposition  of  calcite  HI  or  formation  of  a  phase  not  yet 
identified.  The  final  volumes  of  samples  released  from  the  calcite  m  region  to  P  = 

0  are  large  and  suggest  that  CO2  gas  is  being  evolved  at  very  low  pressures.  There 

is  no  evidence  for  this  behavior  at  P  =  1.41  GPa,  so  that  the  decomposition 
probably  occurs  at  P  <  1  GPa. 

In  contrast  with  release  from  calcite  HI,  the  liquid  phase  seems  to  be 
kinetically  inhibited  from  freezing  down  to  release  state  pressures  as  low  as  4  GPa 
and  follows  the  release  isentropes  described  by  (16).  We  suspect  that  this  phase 
also  undergoes  decomposition  with  evolution  of  CO2  at  P  <  1  GPa. 

As  with  Si02,  we  model  release  from  mixed-phase  regions  of  the  Hugoniot 

as  release  from  the  metastable  Hugoniots  of  the  two  end-member  phases  at  the  same 
pressure.  For  this  purpose,  we  assume  that  release  from  the  calcite  I  metastable 
Hugoniot  is  the  simple  release  isentrope  described  by  (16).  Release  adiabats 
described  by  this  model  are  shown  in  Figure  7. 
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III.  DISCUSSION 


Now  that  a  model  for  the  total  shock-release  cycle  has  been  developed  we 
can  apply  it  to  calculate  the  fraction  f  of  the  shock  internal  energy  that  is 
irreversibly  deposited  in  the  material  for  a  given  shock  pressure  using  the  equation 

f=EH-Ef  (27) 

Eh 

where  EH  and  the  quantity  (Eh-E,)  are  given  by  (4)  and  (17),  respectively.  We  have 
assumed  here  that  the  effects  of  the  elastic  precursor  at  low  pressures  can  be 
ignored.  The  resulting  estimates  for  f  are  presented  in  Figure  8.  For  Coconino 
Sandstone,  f  approaches  1  at  low  shock  pressures  because  of  pore  collapse,  but 
rapidly  drops  to  -0.44  at  18  GPa,  rising  rapidly  to  -0.65  at  60  GPa  and  varying 
slowly  thereafter,  decreasing  to  -0.59  at  200  GPa.  In  the  case  of  the  limestones, 
we  have  ignored  the  rapid  expansion  due  to  C02  gas  evolution  at  extremely  low 
pressures,  so  the  values  of  f  given  are  upper  limits  in  the  case  of  release  to  P=0. 

At  low  pressures,  as  with  sandstone,  f  approaches  1  due  to  pore  collapse,  but 
drops  rapidly,  the  rate  of  the  drop  being  dependent  on  the  initial  density.  Because 
of  differences  in  the  phases  on  the  Hugoniot  (aragonite  in  Bedford  limestone  versus 
calcite  HI  in  Solenhofen  Limestone)  and  the  difference  in  poo,  the  P-f  curves 
behave  very  differently  in  detail  below  120  GPa.  Above  120  GPa,  both  curves  are 
very  similar,  but  with  slightly  greater  energy  deposition  (higher  f )  in  the  more 
porous  Bedford  limestone,  as  should  be  expected. 

There  are  some  fundamental  observations  that  should  be  made  about  the 
curves  in  Figure  8.  The  detailed  structure  of  the  curves  should  be  taken  with  a  fan- 
amount  of  skepticism,  since  the  models  for  the  phase  transition  kinetics  and  the 
phase  boundaries  are  based  on  sparse  data.  However,  the  relative  positions  of  the 
curves  is  a  property  only  of  the  equations  of  state  thus  should  be  a  robust  result. 
Since  a  higher  value  off  is  positively  correlated  with  the  irreversible  energy 
deposition  and  thus  the  rate  of  decay  of  the  shock  wave,  rocks  characterized  by 
high  f  values  will  tend  to  decouple  explosions  from  far  field  observables  more 
efficiently  than  would  rocks  with  low  f  values.  In  particular,  the  conclusion  from 
the  curves  presented  in  Figure  8  is  that  porous  silicate  rocks  such  as  sandstones  are 
much  better  suited  to  decoupling  that  are  carbonate  rocks.  This  is  aside  from  other 
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Figure  8.  Fraction  f  of  Rankine-Hugoniot  internal  energy  which  is  irreversibly  deposited  as  heat  in  Coconino  Sandstone  and 
Solenhofen  and  Bedford  limestones,  shown  as  a  funtion  of  peak  shock  pressure. 


issues,  such  as  the  need  to  contain  gas  evolved  from  decomposition  of  CaCC>3  as 
release  pressures  approach  zero. 

An  important  effect  which  has  been  ignored  here  is  the  effect  of  the  elastic 
precursor  wave  that  precedes  the  plastic  shock  wave  at  low  pressures.  The  actual 
values  of  Eh  for  our  limestone  experiments  which  showed  elastic  precursors  and 
intermediate  waves  were  -35%  and  -65%  for  Bedford  and  Solenhofen  limestones, 
respectively,  of  the  value  that  is  calculated  for  the  same  final  particle  velocity  in  a 
purely  plastic  shock  such  as  we  have  used  in  the  present  calculations.  For  the  data 
from  Ahrens  and  Gregson  [1964]  for  Coconino  Sandstone,  the  actual  internal 
energy  is  -55%  of  the  calculated  plastic  shock  energy.  Also,  material  strength  in 
mitigating  pore  collapse  at  low  pressures,  so  that  the  energy  deposition  is  much 
lower.  Since  the  release  behavior  should  be  the  same,  however,  the  result  is  that  at 
very  low  pressures,  even  in  porous  materials,  f -0.  There  is  ample  experimental 
evidence  that  significant  energy  deposition  is  occurring  at  low  shock  pressures, 
however,  so  that  processes  such  as  fracture  formation  must  be  important  to  the 
overall  question  of  coupling  of  the  near-source  behavior  to  far-field  observables. 
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IV.  CONCLUSIONS 


We  have  performed  a  series  of  shock  Hugoniot  and  release  adiabat 
measurements  on  sandstone  and  limestone  and  combined  the  results  with  previous 
data  to  develop  models  for  the  shock  and  release  behavior  of  these  rocks.  The 
single-phase  portions  of  the  sandstone  Hugoniot  curve  can  be  calculated  from  the 
equations  of  state  of  quartz,  stishovite,  and  a  high  coordination  number  liquid.  The 
properties  of  quartz  and  stishovite  have  been  extensively  studied  in  the  past  We 
present  a  new  fit  of  available  data  for  the  dense  liquid  phase  to  get  poo  =  4.150 
Mg/m3,  Kso  =  313.5  GPa,  K’  =  1.8484,  y  =  1.4(V/V0)1/2,  and  E*  =  2.4  MJ/kg.  In 
the  case  of  Coconino  Sandstone,  the  presence  of  a  few  percent  of  non-quartz 
constituents  requires  that  the  STP  densities  of  the  relevant  phases  be  reduced  by 
2.4%  to  fit  the  data.  The  intervening  mixed-phases  region  can  be  modeled  as  a 
simple  mixture  of  these  two  relevant  end  member  phases  under  the  assumption  that 
formation  of  the  higher  pressure  phase  is  kinetically  inhibited.  The  fraction  f  of  the 
higher  pressure  phase  is  described  by  f  =  Ae  E,/RT.  For  formation  of  stishovite 
from  quartz,  A  =  7.286  and  Ea(kJ  /  mol)  =  89.36  -  71.97(p0/p00)  where  poo  and 

p0  are  the  actual  and  theoretical  nonporous  densities  of  the  initial  rock  sample.  For 
formation  of  the  dense  liquid  from  stishovite,  A  =  6.245  x  10^  and  Ea  (kJ/mol)  = 
404  +  0.312poo.  No  formation  of  the  higher  pressure  phase  occurs  until  the  phase 
boundary  is  reached.  As  crude  approximations  we  describe  the  quartz-stishovite 
phase  boundary  by  P  =  12  GPa  and  the  stishovite-dense  liquid  phase  boundary  by 
T  =  4348  -  12.121P  for  T  in  K  and  P  in  GPa. 

Release  adiabats  from  the  single  phase  regions  of  the  sandstone  hugoniot 
are  modeled  by  the  isentropes  of  the  relevant  phases,  as  described  by  the  equation 
of  state.  At  ~5  GPa,  stishovite  begins  to  form  a  diaplectic  glass.  We  model  the 
release  adiabat  below  this  pressure  as  a  straight  line  in  P-V  space,  from  the 
stishovite  release  isentrope  at  5  GPa  to  the  STP  diaplectic  glass  density  of  2.27 
Mg/m3.  The  dense  liquid  undergoes  a  similar  change  beginning  at  ~5  GPa  with 
formation  of  normal  silica  glass  with  an  STP  density  of  2.204  Mg/m3.  The  quartz 
release  adiabat  follows  the  quartz  isentrope  to  P=0.  As  with  the  Hugoniot,  we 
model  the  release  from  the  mixed-phase  regions  of  the  Hugoniot  as  a  mixture  of  the 
relevant  phases,  with  the  proportions  frozen  at  the  Hugoniot  values. 

The  limestone  Hugoniot  curve  is  complicated  by  several  metastable  phases 
and  multiple-wave  shock  structure  at  low  pressures.  Single  phase  portions  of  the 
Hugoniots  of  Solenhofen  and  Bedford  limestones  are  modeled  as  calcite  I  at  low 
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pressures  and  a  dense  liquid  at  high  pressures.  We  have  developed  an  equation  of 
state  for  the  dense  liquid  phase,  with  p0  =  3.377  Mg/m3,  Kso  =  139  GPa,  K'  =  3.4, 
y  =  2.5(V/V0),  and  E*  =  2.142  MJ/kg.  At  intermediate  pressures,  dense  limestones 
such  as  Solenhofen  seem  to  form  calcite  HI,  while  the  more  porous  Bedford 
limestone  apparently  forms  aragonite.  The  mixed  phase  regions  of  the  Hugoniot 
are  modeled  in  a  manner  analogous  to  that  used  for  the  sandstone  Hugoniot.  The 
parameters  describing  the  fraction  f  of  the  higher  pressure  phase  in  the  mixed  phase 
regions  are  A  =  5.848  and  Ea  =  62.713-20.651p00  for  formation  of  calcite  ffl  or 
aragonite  from  calcite  I  and  A  =  2.919  and  Ea  =  382.453-  130.448Poo  for  formation 
of  the  dense  liquid  from  either  aragonite  or  calcite  III,  for  Ea  in  kJ/mol  and  poo  hi 
Mg/m3.  The  phase  boundaries  in  our  model  are  approximated  by  T^  =  164.7P- 
3-6,  Tun  =  78.00-1.3,  =  745+48P,  and  Tin.Liq  =  401+25.8P. 

As  with  sandstone,  we  model  the  release  adiabats  with  the  isentropes  of  the 
relevant  phases,  mixing  them  at  the  appropriate  ratios  for  release  from  the  mixed- 
phase  portions  of  the  Hugoniot.  There  is  no  evidence  of  formation  of  a  glass  at 
very  low  pressures,  but  release  of  phases  other  than  calcite  I  is  apparently 
accompanied  by  decomposition  with  evolution  of  C02  gas  at  pressures  below  1 
GPa.  We  have  ignored  this  very  low  pressure  effect  in  our  calculations. 

Use  of  these  results  allows  us  to  estimate  the  fraction  f  of  shock  internal 
energy  irreversibly  deposited  in  these  rocks.  This  energy  in  no  longer  available  to 
drive  the  shock  wave,  so  the  shock  wave  decays  more  rapidly  than  expected  from 
geometric  effects  alone.  We  find  that  f  is  a  factor  of  ~2  smaller  in  limestones  than 
in  sandstone,  so  that  decoupling  of  the  source  from  far-field  seismic  waves  is  more 
efficient  in  sandstones.  This  analysis  ignored  the  effects  of  material  strength  at  low 
shock  pressures,  although  such  effects  are  important.  The  value  of  f  may  be 
reduced  considerably  in  the  presence  of  an  elastic  precursor  wave  and  of 
intermediate  waves.  Experimental  evidence  indicates  that  shock  wave  decay  still 
occurs  even  in  the  supposedly  purely  elastic  shock  regime,  showing  that  effects 
such  as  fracture  formation  are  very  important  to  the  overall  picture  of  explosive 
energy  coupling. 
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